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1. Introduction 

RHIC collisions are commonly described in terms of two themes: hydrodynamic (hydro) evo- 
lution of a thermalized bulk medium and energy loss of energetic partons (hard probes) in that 
medium. Hydro is thought to dominate p t spectra below 2 GeV/c, parton fragmentation is expected 
above 5 GeV/c, and "quark coalescence" is thought to dominate the intermediate p t interval. 

Recent studies of spectrum and correlation structure have revealed interesting new aspects of 
RHIC collisions. Number and p t angular correlations in the final state contain minijet structures 
(minimum-bias parton fragmentation) [fil @, 0, EL BL E| fj\\. Two-component analysis of p-p and 
Au-Au spectra reveals a corresponding hard component interpreted as a minimum-bias fragment 
distribution, suggesting that jet phenomena extend down to 0.1 GeV/c hadron momentum 



Minijets (well described in p-p collisions by PYTHIA/HIJING [10]) are observed to dom- 
inate the transverse dynamics of nuclear collisions at energies above ^/snn ~ 15 GeV. The term 
"minijets" can be applied collectively to hadron fragments from the minimum-bias scattered-parton 
spectrum averaged over a given A-A or N-N event ensemble. Minijets provide unbiased access to 
fragment distribution structure down to a small cutoff energy for scattered partons (those partons 
fragmenting to charged hadrons) and to the smallest detectable fragment momenta (~ 0.1 GeV/c). 

In this analysis minijets are studied in the form of p r spectrum hard components isolated via 
the two-component spectrum model. Measured hard components are compared with calculated 
fragment distributions obtained by folding parton spectra with fragmentation-function ensembles. 
Parton spectrum parameters and modifications to fragmentation functions in more-central Au-Au 



collisions are inferred [11]. The goal is a comprehensive QCD description of all nuclear collisions. 



2. Two-component spectrum model 

The two-component model of p-p spectra [Q] is the starting point for the fragmentation analy- 
sis described here. The two-component (soft+hard) model was first obtained from a Taylor-series 
expansion on observed event multiplicity n c /, (< corrected n c h) of spectra for several multiplicity 
classes. The soft component was subsequently interpreted as longitudinal projectile-nucleon frag- 
mentation, the hard component as transverse scattered-parton fragmentation. The two-component 
model applies to two-particle correlations on (y t ,y>t) as well as their ID projections onto p t or y t . 

The two-component spectrum model for p-p collisions with corrected soft and hard multiplic- 
ities n s + nh = n c h is 

1 ldnch ^-S (y t ) + n J^H (y t ), (2.1) 



n s {h ch )yt dy t n s {n ch ) 

where soft component So(y t ) is the Taylor series "constant," and hard component Ho(y t ) is the 
coefficient of the term linear in n c h, both normalized to unit integral. So(y t ) is a Levy distribution on 
m t , Ho(y t ) is a Gaussian plus QCD power-law tail on transverse rapidity y t = ln{ (m t +p t )/mo}. To 
compare with A-A spectra we define S pp = (1 /y t )dn s /dy t with reference model h s Sq and similarly 
for H pp <-> nh Ho. The two-term Taylor series exhausts all significant p-p spectrum structure. 

Fig. [j] (first panel) shows spectra for ten multiplicity classes from 200 GeV non-single diffrac- 
tive (NSD) p-p collisions [ph. The asymptotic limit for h c h — » (dash-dotted curve) is Sq. The 
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spectra are normalized by soft-component multiplicity n s . Fig. [I] (second panel) shows the two- 



component algebraic model Eq. (gjj) with unit-normal model functions Sq and Ho defined in | 
Hard-component coefficient n n jn s scales as an cn . Factor a = 0.01 is the average value for most 
n cn classes. The spectrum data in the first panel are described to the statistical limits. 
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Figure 1: First: y t spectra from ^/snn — 200 GeV p-p collisions for ten multiplicities, Second: Correspond- 
ing two-component model, Third: Corresponding hard components, Fourth: Hard components normalized 
to NSD p-p collisions. 



Figure [I] (third panel) shows hard components H pp /n s for ten multiplicity classes obtained by 
subtracting fixed soft component Sq from the ten NSD p-p spectra normalized to n s . The shape is 
Gaussian independent of multiplicity [§]. Fig. [j] (fourth panel) shows hard components H pp from 
the third panel scaled by factors nh(\.25)/tih{n c h) to obtain the mean hard component for NSD 
p-p collisions. The dash-dotted curve is 0.02 # [0.02 ~ (a = 0.007) (h ch = 1.25) (n s = 2.5) [§]]. 
The exponential tail represents the QCD power law <=c p t " eCD . The spectrum hard component 
is interpreted as a minimum-bias fragment distribution dominated by "minijets" — jets from those 
partons (gluons) with at least the minimum energy required to produce charge-neutral combinations 
of charged hadrons. Equivalent structure appears in two-particle correlations on (y t ,y>t) 

The corresponding two-component model for per-participant-pair A-A spectra is 

2 1 dn c h / \ , 

" ■■ Snn^ + vHaa^v) (2.2) 



ripanyt dy t 

= SNN(yt) + VrAA(y t ;v)H NN (y t ), 

where Snn (~ S pp ) is the soft component and Haa is the A-A hard component (with reference 
Hnn ~ H pp ) integrating respectively to multiplicities n s and rih in one unit of pseudorapidity rj [|| 
^]. Ratio yaa = Haa/ Hnn is an alternative to nuclear modification factor Raa- Centrality measure 
V = 2nt,inary/n participant estimates the Glauber-model mean nucleon path length. We are interested 
in the evolution of hard component (fragment distribution) Haa or ratio taa with A-A centrality. 



3. Fragmentation functions 

e + -e~ (e-e) fragmentation functions (FFs) have been parametrized accurately over the full 
kinematic region relevant to nuclear collisions, e-e light-quark and gluon fragmentation functions 
D X x(x,Q 2 ) *-* D xx {y,y max ) (xx is the FF context: e-e, p-p, A-A) are accurately described above 
energy scale (dijet energy) Q ~ 10 GeV by a two-parameter beta distribution ft(u;p, q) on normal- 
ized rapidity u Jl2[]. Fragment rapidity for unidentified hadrons is y = \n[(E + p)/m n ], and parton 
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rapidity y max = \n{Q/m n ). Parameters (p,q) vary slowly and linearly with y max above 10 GeV and 
can be extrapolated down to Q ~ 4 GeV based on dijet multiplicity data. 

Fig. || (first panel) shows measured FFs for three energy scales from HERA/LEP [13, 
The 2 in the axis label indicates dijet n c h densities. The vertical lines at right denote y max values. 
The curves are determined by the j3(p,q) parametrization with y m j„ ~ 0.35 (p t ~ 0.05 GeV/c, left 
vertical line) and describe data to their error limits over the entire fragment momentum range. 
Fig. U (second panel) shows the FF ensemble (inclusive light quarks fragment to inclusive hadrons) 
vs energy scale Q as a surface plot [12]. The dashed curve is the locus of modes — the maximum 
points of the FFs. Between the dash-dotted lines the system is determined by FF data. Between the 
dash-dotted and dotted lines the parametrization is constrained only by dijet multiplicities. 




10 , 10 
Q = (GeV) 

Figure 2: First: Fragmentation functions (FFs) from e + -e~ collisions for three energies with j3 -distribution 
parametrizations (solid curves), Second: Full e + -e~ FF parametrization on parton rapidity y max , Third: FFs 
from p-p collisions for several dijet energies, Fourth: Full p-p FF parameterization on parton rapidity. 

Figure || (third panel) shows FF data from p-p collisions at FNAL (samples from the full data 
set) [|TJ]. The solid curves guide the eye. There is a significant systematic difference between 
p-p and e-e FFs. The dotted line represents the lower limit for e-e FFs. The systematic gap for 
all parton energies is apparent — y m ,„ for p-p collisions is ~ 1.5 (0.3 GeV/c) instead of 0.35 (0.05 
GeV/c). The CDF FFs also reveal a systematic amplitude saturation or suppression at larger parton 
energies compared to LEP systematics. The curve labeled MB is the hard-component reference 
from NSD p-p collisions [JSp . Fig. ^(fourth panel) shows a surface plot of the p-p FF ensemble [p"T|]. 
The surface represents the e-e FF parametrization modified by introducing cutoff factor 



gem (y) = tanh{ (y - y ) (E, y } y>y , 



(3.1) 



with y ~ £ y ~ 1.5 determined by the CDF FF data [[To]]. The modified FFs have not been rescaled 
to recover the initial e-e parton energy. The cutoff function thus represents real fragment and energy 
loss from p-p relative to e-e FFs. The difference implies that FFs are not universal. 

Figure || (first panel) shows parametrized beta FFs for five e-e energy scales. The Q = 6 GeV 
scale is associated with minijets as explained below. Such curves provide a complete description 
of e-e FFs at energy scales relevant to nuclear collisions. Fig. ^ (second panel) shows light-quark 
dijet multiplicity systematics from the same beta parametrization. The solid points correspond to 
the FFs in the first panel. The open circles represent multiplicities from medium modification of 
those FFs in central Au-Au collisions at 200 GeV, as described in Sec. ^. The "in-medium" shift of 
FFs to smaller fragment momenta requires more fragments to satisfy parton-energy conservation. 
The systematics of quark and gluon jets coincide for energy scales Q = 2Ej et < 8 GeV (y max < 4). 
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Figure 3: First: Parametrized e + -e~ FFs for five dijet energies, Second: Corresponding dijet multiplicities 
for in-vacuum (solid points) and in-medium (open points) FFs, Third: Parametrized p-p FFs for five dijet 
energies compared to CDF data (points) [n5[j, Fourth: Corresponding dijet multiplicities for p-p FFs (solid 
points) and published values (open points [O, |l8||). 



Figure || (third panel) shows e-e beta FFs for five parton energies [ |r2| ] modified by the g cut 
factor to describe p-p FFs. The deviation from e-e FFs is indicated by the two dotted lines JTl|]. 
The CDF data (points) are from [15]. Fig. || (fourth panel) shows multiplicity systematics (solid 
points) for p-p (i.e., modified e-e) FFs from the parametrization. The solid curve represents un- 
modified e-e FFs as a reference. There is substantial reduction of p-p FF multiplicities due to the 
cutoff. Also plotted are CDF FF multiplicities from reconstructed jets (open triangles [|n|] and 
open circles Jl8|]). 

Comparison of Fig. || second and fourth panels reveals that dijet multiplicities (and charged- 
particle energy integrals) are strongly suppressed in p-p collisions compared to equivalent FFs in 
e-e collisions, p-p jet multiplicities are reduced by 30-70%. FFs are apparently "modified" in p-p 
collisions as well as A-A collisions. At Q = 6 GeV (minijets) there is a three-fold dijet multiplicity 
reduction for p-p relative to e-e collisions. 



4. Parton spectrum model 

A model for the parton p t spectrum resulting from minimum-bias scattering into an r/ accep- 
tance near projectile mid-rapidity can be parametrized as 



1 do d , 



ijet 



: fcut (Pt, 



ijet 



fcut (ymax) Ay max exp{-(n QC D-2)y ma x}, (4.1) 



Pt dp t p t QCD dy max 

which defines QCD exponent uqcd, with y max = ln(2 p t /m n ). The cutoff factor 

fcutiymax) = {tanh[0 

max ycut )/&*] + 1}/2 



(4.2) 



represents in this analysis the minimum parton momentum which leads to detectable charged 
hadrons as neutral pairs (i.e., local charge ordering). Parton spectrum and cutoff parameters are 
determined via FD comparisons with p-p and Au-Au spectrum hard components. 

Fig. ^| (semilog and linear formats) shows the parton spectrum (solid curve) inferred from a 
p-p spectrum hard component [O]. y cut and A ymax are well-defined by the p-p hard component, 
and uqcd is defined by Au-Au spectrum hard components extending to larger y t . The dotted curve 



in the first panel is an ab-initio pQCD calculation [19]. The linear plot (second panel) indicates 
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Figure 4: First: Parton spectra inferred from this analysis for p-p collisions (solid curve) and central Au-Au 
collisions (dash-dotted curve) compared to an ab-initio pQCD theory result (bold dotted curve [|l9]|), Second: 
Parton spectra from this analysis in a linear plot, Third: Parton spectrum from reconstructed jets (UA1, solid 
points [^0|) compared to theory (dashed curve [|2j]) and this analysis (solid curve, note factor 3). 



the narrowness of the spectrum, with effective mean energy near 3 GeV (minijets). Fig. |] (third 
panel) compares the spectrum defined in this analysis (solid curve, and note the factor 3) with 200 



GeV UA1 jet cross-section data obtained by event- wise jet reconstruction [2C]. The UA1 spectrum 



integral is 4 mb |2C]. The spectrum from this analysis integrates to 2.5±0.6 mb with well-defined 
cutoff ~ 3 GeV which agrees well withpQCD theory (e.g., [21]). The KLL parametrization 600/ p] 



mb/(GeV/c) (dashed line) integrates to 2.2 mb above 3 GeV/c [22]. 



5. Fragment distributions from a QCD folding integral 

The folding integral used to obtain fragment distributions (FDs) in this analysis is 
d 2 n h ^ 6(5t],At}) n , ^dOdtjet 



, , , dymaxD^y^ymax) — , (5.1) 

dydr] o NSD Ar] Jo dy max 

where Acc^ymox) 1S the dijet FF ensemble from a source collision system (xx = e-e, p-p, A-A, 



in-medium or in- vacuum), and do c uj et /dy max is the minimum-bias parton spectrum Ql 1Q. Hadron 
spectrum hard component d 2 nh/dydr\ as defined represents the fragment yield from scattered par- 
ton pairs into one unit of rj. Efficiency factor e ~ 0.5 (for a single dijet and one unit of Tj) includes 
the probability that the second jet also falls within v\ acceptance 8f] and accounts for losses from 
jets near the acceptance boundary. At] ~ 5 is the effective 471 v\ interval for scattered partons. Onsd 
(~ 36 mb for y/siw = 200 GeV) is the cross section for NSD p-p collisions. 



Fig. H| (first panel) shows the integrand D ee (y, y max ) d ^ >! " of the folding integral in Eq. ( JsTT] ) in 



corporating unmodified FFs from e-e collisions with lower bound at y m ;„ ~ 0.35 (p t ~ 0.05 GeV/c) 
(dotted line). The plot z axis is logarithmic to show structure over the entire distribution support. 
Fig. H (second panel) shows the corresponding FD (solid curve). The parton spectrum parameters 
determined from the p-p hard component are retained. The solid curve is the "correct answer" for 
an FD describing inclusive hadrons from inclusive partons produced by free parton scattering from 
p-p collisions, which is not observed in real nuclear collisions. The dash-dotted curve represents 
the hard-component model inferred from p-p collisions [^]. The FD from e-e FFs lies well above 
the measured p-p hard component for hadron p t < 2 GeV/c (y t < 3.3), and the mode is shifted 
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down to ~ 0.5 GeV/c. The "correct" e-e FD strongly disagrees with the relevant part of the p-p p t 
spectrum — the hard component. Despite the strong disagreement the e-e FD is the correct reference 
for nuclear collisions, as demonstrated below. 




Figure 5: First: pQCD folding-integral argument for e + -e~ FFs, Second: Corresponding fragment distribu- 
tion (solid curve) and p-p hard-component reference (dash-dotted curve), Third: Folding-integral argument 
for p-p FFs, Fourth: Corresponding fragment distribution (solid curve) compared to p-p hard-component 
data (points). Dotted curves correspond to ±10% change in parton spectrum cutoff energy about 3 GeV. 

Fig. R (third panel) shows a surface plot of integrand D pp (y,y max ) d ^ d, ' e \ incorporating e-e 
FFs based on the LEP parametrization but modified by the FF cutoff function inferred from p- 
p collisions. The main difference from e-e FFs is that the lower bound of p-p FFs is raised to 
ymin ~ 1-5 (p t ~ 0.3 GeV/c from 0.05 GeV/c). Fig. [5] (fourth panel) shows the corresponding FD 
Hnn-vcic (integration of the third panel over y max ) as the solid curve. The mode of the FD is ~ 1 
GeV/c. The dash-dotted curve is a Gaussian-plus-tail model function, and the solid points are hard- 
component data from p-p collisions [|8|]. That comparison determines parton spectrum parameters 
ycut = 3.75 {E cut ~ 3 GeV), A ynm and exponent uqcd = 7.5. The p-p data are well-described by 
the pQCD folding integral. This procedure establishes an absolute quantitative relationship among 
parametrized parton spectrum, measured FFs and measured spectrum hard components over all p t , 
not just a restricted interval (e.g., above 2 GeV/c). 



6. Parton "energy loss" and medium-modified FDs 

The hypothesis of parton energy loss in a thermalized bulk medium is of central importance at 
RHIC. In some models the medium is opaque to most hard-scattered partons - only a small fraction 
emerge as correlated fragments. But minijet systematics suggest no parton loss to thermalization. 
In this section I adopt a pQCD-inspired minimal model of FF modification (Borghini-Wiedemann 
or BW) [E3], with no loss of parton energy to a medium or scattered partons to thermalization. 



Figure g (first panel) illustrates the BW model of FF modification (cf. Fig. 1 of []23|]). In- 
vacuum e-e FFs for Q = 14 and 200 GeV from the beta parametrization are shown as dashed and 
solid curves respectively [12]. Whereas the BW model was expressed on £ p FFs are plotted here 



on fragment rapidity y. The relation is E, p = \n{pj et /p) = ln(2 pjet/m^) — \n{2p/m n ) ~ y max — y, 
with energy scale Q = 2pj et . The practical consequence of the BW "energy-loss" mechanism is 
a momentum-conserving rescaling of FFs on x p , with ^ p = ln(l/x p ). Small density reductions 
at larger fragment momenta (smaller £ p ) are compensated by much larger increases at smaller 
momenta. The largest changes (central Au-Au) correspond to an inferred 25% leading -parton frac- 
tional "energy loss." We model the BW modification simply by changing parameter q in f}(u;p,q) 
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by Aq ~ 1, which accurately reproduces the BW result. The modified FFs are the dash-dotted and 
dotted curves JTT|]. Fig. || (second panel) shows the modified e-e FF ensemble with FF modes 
shifted to smaller fragment rapidities. No energy is lost from FFs in this model. 




2 4 6 
y = ln[(E + p)/m ] 



Figure 6: First: e + -e FFs for two energies unmodified (solid and dashed curves) and modified according 
to a rescaling procedure [0 (dash-dotted and dotted curves) to emulate parton "energy loss," Second: e + - 



e~ FF ensemble modified according to p3[ , Third: Medium-modified FD from e + -e~ FFs (solid curve) 
compared to in-vacuum e + -e~ FD (dotted curve) Fourth: Medium-modified FD from p-p FFs (solid curve) 
compared to in-vacuum FD (dotted curve). 

Figure |6| (third panel) shows H ee ^ me( j (solid curve), the FD obtained by inserting e-e in-medium 
FFs from the second panel into Eq. ( |Q| ) and integrating over parton rapidity y max . The dotted curve 
is the H ee - vac reference from in-vacuum e-e FFs. The dash-dotted curve is again the Gaussian-plus- 
tail p-p hard component Hqq reference. The mode of H ee - me d is ~ 0.3 GeV/c. Fig. ^ (fourth panel) 
shows results for p-p FFs. Major differences between p-p and e-e FDs appear below p t ~ 2 GeV/c 
(y t ~ 3.3). Conventional comparisons with theory (e.g., data vs NLO FDs) typically do not extend 



below 2 GeV/c [24]. The large difference between the two systems below 2 GeV/c reveals that the 
small- p t region, conventionally assigned to hydro phenomena, may be of central importance for 
understanding fragmentation evolution in A- A collisions. 



7. Fragment evolution with centrality in Au-Au collisions 

We have established a system to combine measured FFs and a parametrized pQCD parton 
spectrum to produce calculated fragment distributions FD XX for comparison with measured spec- 
trum hard components H xx . Conventional comparisons employ a ratio measure. Two questions 
emerge: what is the validity of the ratio definition, and what should be the reference for such a 
ratio. The conventional spectrum ratio at RHIC is Raa, defined in the first line of 

S NN (y t ) + vHAA(y t ;v) 



R 



AA 



X 



- + T r AA 



+HNN(yt) 

at y t = 2. 



(7.1) 



In that definition the terms in numerator and denominator are normalized per participant pair 
n part /2, so the prefactor is 1/v rather than l/nbj nary . Fig. [7] (first panel) illustrates problems with 
that measure. Hard-component evolution with centrality, the main object of this analysis, is de- 



scribed by ratio r&A = Haa/FInn- The second line of Eq. (7.1) gives the limiting value of Raa near 
y t ~ 2 where the Hnn/S^n ratio is typically ~ 1/170. r^A is thus suppressed by a large factor 
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in just the interval where fragmentation details are most important. The p-p data (dots) illustrate 
suppression of even statistical fluctuations. All information is lost. 




Figure 7: First: Conventional spectrum-ratio measure Raa, illustrating strong suppression of specttum 
information below 4 Ge V/c (y, — 4), Second: Hard-component ratio r^ illustrating restoration of suppressed 
sttucture at small y t , Third: Comparison of calculated FD ratios to measured for central Au-Au collisions, 
Fourth: Comparison of novel FD ratio r e N to measured r^A for central Au-Au collisions. 

Figure ^ (second panel) shows ratio yaa based on hard-component reference H^n set equal to 
Gaussian model Hqq = rif,Ho from Evolution of suppression and enhancement is dramatically 
more accessible. The p-p data and the most peripheral Au-Au data agree with the N-N reference 
( r AA = 1) above y t = 2.5 but deviate significantly from H GG below that point. For the Au-Au 
collisions in this figure v = 2nbi n /n part values for five centralities are 1.93, 2.83, 3.92, 4.87, 5.5, 
where v ~ 1.25 is N-N collisions and V ~ 6 is b = Au-Au collisions From v = 1.98 to v = 
2.83 there is a dramatic change in the hard component. At the transition point v ~ 2.5 n part = 40 
(out of 382) and n bin = 50 (out of 1 136). 

Figure [7] (third panel) shows calculated FD ratios r x 



FD xx - med /FD xx „ vac with xx = e-e 



(dash-dotted curve, e-e FFs) or N-N (dashed curve, p-p FFs) [|11|]. The solid curve is the measured 
yaa from central (0-12%) Au-Au collisions at 200 GeV Aq ~ 1.15 for H ee - me d and //Aw-merf 
(in-medium FFs) was adjusted to obtain the correct large-y ? suppression for 0-12% central Au-Au. 
The reference for r^A is hard-component model function Hqg- The dotted curve is a reference ratio 
obtained by shifting Hqg on y t by Ay t ~ —0.26 (negative boost) [|p. The simple negative-boost 
model does not describe the Au-Au data. But the e-e and N-N ratios also do not describe the data. 

Figure [7| (fourth panel) introduces a novel concept. Instead of comparing the calculated in- 
medium FD for N-N collisions averaged within A-A collisions with the in-vacuum FD for isolated 
N-N collisions, or similarly comparing e-e with e-e as in the third panel, the in-medium FD for e-e 
is compared with the in-vacuum FD for N-N by defining ratio 

FDee— med 



r e N 



(7.2) 



FDpjN-vac 

Calculated r e N describes the measured taa well over the entire fragment momentum range. We 
conclude that FDnn-vuc is not the correct reference. The proper in-vacuum reference for all sys- 
tems is an FD from e-e FFs, not p-p FFs. We define FD ratios r xx = FD xx ^ yyy /FD ee ^ vac with xx 
= ee, NN, AA and yyy = med or vac to be compared with equivalent spectrum hard components 

H xx —yyy ■ 

Figure || (first panel) shows ratios redefined in terms of the ee-vac reference: H pp (p-p data 
- points), Haa (peripheral Au-Au data - solid curve) and calculated H ee _ me d (dash-dotted curve) 
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and HNN- vac (dashed curve) all divided by reference H ee - vac . The strong suppression of p-p and 
peripheral Au-Au data apparent at smaller y t results from the cutoff of p-p FFs noted above. The 
comparison is linear rather than logarithmic, as in Fig. [7|, and is thus more differential. 




Figure 8: First: FD ratios relative to an ee-vacuum reference for Au-Au collisions below the sharp transi- 
tion, Second: FD ratios relative to an ee-vacuum reference for Au-Au collisions above the sharp transition 
revealing major changes in FD structure, Third: Hard-component evolution in central Au-Au collisions vs 
centrality. Large increases in fragment yield at smaller y t (p t < 2 GeV/c) accompany suppression at large y t . 



Figure [8] (second panel) shows measured H^jH ee ^ vac for more-central Au-Au collisions (solid 
curves) above a transition point on centrality at v ~ 2.5. The main difference is partial restoration 
of the suppressed region at smaller y t and suppression at larger y t . The latter has been the ma- 



jor observation at RHIC for jet-related modification (high-p f suppression, "jet quenching" [25]). 
Apparent from this analysis is the accompanying very large increase in fragment yield below 2 
GeV/c, still strongly correlated with the parent parton [J7|]. Also notable is the substantial gap be- 
tween the peripheral data and the four more-central spectra [§]. Changes in fragmentation depend 
very strongly on centrality near the transition point. It is remarkable that the trend at 10 GeV/c 
corresponds closely to the trend at 0.5 GeV/c. Calculated FD ratio r ee (dash-dotted curve) cor- 
responds to a parton spectrum cutoff shifted down to 2.7 GeV from 3 GeV for p-p collisions, as 
shown in Fig. ^| (first and second panels). The shift may result from an increased hadron density of 



states [11]. 



Figure || (third panel) shows spectrum hard components Haa (solid curves) for five centralities 
from 200 GeV Au-Au collisions [^]. The hard components of y t spectra scale proportional to 
ribinary, as expected for parton scattering and fragmentation in A-A collisions (jets). The points are 
hard-component data from 200 GeV NSD p-p collisions [j|]. The dash-dotted curve is the standard 
Gaussian+tail model function Hqg- Calculated FDs are also shown. The dashed curve is H^N-vac, 
and the upper dotted curve is fi ee — mec i with Aq = 1.15, which nominally corresponds to the most- 
central Au-Au curve (0-12%). The parton spectrum cutoff for H ee - mec i has been reduced from 3 
GeV {y max = 3.75) to 2.7 GeV (y max = 3.65) to match the central Au-Au hard component near 
y t = 3. The dotted curves labeled 2 and 5 (Au-Au centralities) are H ee - mec i with cutoff parameters 
yo = E, y reduced to accommodate the data below y t = 2.5. H pp , Haa and ratios based on the e-e 



in- vacuum reference are thus well described by pQCD FD ratio data from 0.3 to 10 GeV/c [11] 
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8. Discussion 

This analysis establishes a quantitative correspondence between calculated pQCD FDs and 
measured spectrum hard components H xx over the entire fragment p t range and parton spectrum. 
We obtain direct access to medium-modified FFs and the underlying parton spectrum. 

In p-p and in peripheral Au-Au collisions below a transition point at v ~ 2.5 the underly- 
ing power-law parton spectrum terminates near 3 GeV. Hard component H pp or Haa is strongly 
suppressed at smaller y t (jet bases excluded from the acceptance) corresponding to p-p FFs. The 
suppression mechanism may be hard-Pomeron (color singlet) exchange in N-N collisions leading 
to color connections different in p-p than in e-e collisions (which produce q-q color dipoles). 

Above the transition point: 1) Measured Haa is strongly enhanced at smaller y t (FF bases 
partially restored) but suppressed at larger y t (so-called "jet quenching"), as observed in [^]. 2) 
Corresponding calculated FDs can be generated by incorporating a "medium-modified" e-e FF 
scenario — simple rescaling of e-e splitting functions — which implies a three-fold increase in jet 
multiplicity compared to p-p FFs. 3) The parton spectrum cutoff is reduced, by up to 10% in 
central Au-Au collisions implying a 50% increase in the jet cross section and minijet production. 

Evolution of Haa corresponds to two-particle correlations on (y t ,yt) [01- Observed spectrum 
hard-component systematics indicate that no partons are "absorbed" or lost to thermalization (no 
"opaque core" is formed). All scattered partons predicted by a pQCD differential cross section 
produce jet-correlated hadrons in the final state. The minimum-bias jet fragment yield in central 
Au-Au collisions fully accounts for the increase of collision multiplicity beyond participant scaling 
(soft component). There is also no indication from correlations, spectrum structure or integrated p t 
that parton spectra extend down to 1 GeV as suggested by saturation-scale arguments [2~6| ], 

9. Summary 

Two-component decomposition of hadron spectra from p-p and Au-Au collisions isolates 
minimum-bias parton fragment distributions as spectrum hard components (// xv ) which can be 
estimated theoretically by folding measured fragmentation functions (FFs) with a pQCD parton 
spectrum to produce calculated fragment distributions (FDs). In this analysis accurate parameteri- 
zations of p-p and e + -e~ FFs for a large range of parton energies are folded with a power-law parton 
spectrum with cutoff to produce calculated FDs which are compared with measured spectrum hard 
components from p-p collisions and from Au-Au collisions for several centralities. 

Comparisons reveal that FFs in p-p collisions are strongly suppressed for smaller fragment 
momenta (jet base suppressed). The suppression is possibly related to hard-Pomeron exchange and 
resulting color-field deviations from q-q. Comparisons further indicate that above a specific Au- 
Au centrality (transition point) there is evolution toward e-e FFs as an asymptotic limit (jet base 
partially restored). FFs are modified consistent with alteration of parton splitting. No partons are 
lost to absorption or thermalization (no "opaque core"), and no significant parton energy is lost 
from integrated FFs. Perturbative QCD describes parton scattering and fragmentation in nuclear 
collisions over a large kinematic domain, and minijets dominate collision dynamics in all cases. 
The most dramatic alteration of parton fragmentation in A-A collisions occurs below p t = 2 GeV/c. 
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